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We study the possibility of switching the types of symmetry breaking bifurcation (SBB) in the 
cyhnder shell waveguide with helical double- well potential along propagation direction. This model 
is described by the one-dimensional nonlinear Schrodinger (NLS) equation. The symmetry- and 
antisymmetry-breakings can be caused by increasing the applied voltage onto the waveguide in the 
self-focusing and -defocusing cases, respectively. In the self-focusing case, the type of SBB can be 
switched from supercritical to subcritical. While in the self-defocusing case, the type of SBB can 
not be switched because only one type of SBB is found. 
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Phase transition is a fundamental topic in many branches of physics. Spontaneous symmetric breaking (SSB), 
which is the transition from the symmetric ground state (GS) to that which does not follow the symmetry of the 
potential, always plays an important role in inducing the transition In nonlinear systems, because additional 
nonlinearity always give rise to the effect of SSB, phase transition and SSB are also important issues and attract great 
O I attention. Among many kinds of nonlinear systems, double-well potential (DWP) or dual-core system is the most 
O , essential model employed to study the phase transition and SSB of the nonlinear states In DWP system, the 

important process relating to the phase transition and SSB is symmetric breaking bifurcation (SBB), which determines 
the process of symmetric states transiting to the asymmetric ones There are two kinds of SBB, subcritical 

' ^ ■ and super critical, which are tantamount to the phase transition of the first and second kind respectively. The former 
^-^1 kind of SBB, subcritical type, features branches of the asymmetric states going at first backward after the bifurcation 
. point then turning forward. While the latter one, the supercritical type, features the asymmetric branches going 
, immediately to the forward direction after the bifurcation point. Different settings with different environments may 

possibly lead to different types of SBB. An interesting problem is the possibility to control the type of the SBB, and 
thus to switch between the respective phase transitions of the two kinds. About this problem, it has been reported 
that the addition of a periodic potential (optical lattice) acting in the unconfined direction changes the character 
of the SBB from sub- to supercritical [7]. Recently, it is also reported that the periodical modulation of the linear- 
coupling constant, which is induced by the electromagnetically induced transparency (EIT) via the double- A system, 
can change the SBB from sub- to supercritical with the increase of the total power of the probe beams p^Tf . 

Very recently, a setting with rotational (alias helical) DWP potential in a cylinder shell waveguide was performed 
[l8|, which reported that the rotating speed of the potential could also induce the symmetry breaking of the non- 
linear modes. It is well known that rotating systems often provide more interesting dynamics than their stationary 
counterparts [l9|-[26|. A natural problem is, whether we can control the type of SBB by means of the rotating speed 
of the potential or not? 

In this work, we build a cylinder waveguide shell with a helical DWP by means of the biased photorefractive crystals 
(BPRCs) (See in Fig. [1]) to study the controllability of type of SBB in the rotating system. The BPRC, which features 
a saturable nonlinearity, is a kind of important nonlinear optical material, especially in creating the spatial optical 
solitons. In the past decades, many kinds of spatial solitons are observed from BPRCs |27l-l3^. As illustrated in 
the Ref. [l8[, helical DWP in the waveguide shell can be created by means of optical-induction techniques. For 
example, using two waves, which have the same amplitude but opposite vortical pumps, ordinary polarizations in 
a photorefractive material, and a very small mismatch of the propagation constant, can form a helical DWP. If the 
width of the shell is thin enough, the problem of this type can be reduced to a quasi one-dimensional (ID) with the 
argument chosen to be the azimuth angle 6>, which require the equation obeying the periodical boundary condition. 
The propagation dynamics of the probe wave, which is launched in the extraordinary polarization of the BPRC, can 
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FIG. 1: Biased photorefractive cylinder waveguide shell, ro is the radius of the circular and U is the applied biased voltage 
onto the crystal. 



be described by the underlying ID scaled nonlinear Schrodinger equation in the rotating reference frame, 



En 



1 



(1) 



where ro, which can be fixed to 1, is the radius of the cylinder shell, uj is the rotating speed of the potential, / is the 
induction field (alias the pumps) and Eq is the nonlinear parameter which is defined by the magnitude and polarity 
of the applied biased voltage U. The periodical boundary condition requires the probe and the pump satisfying 
iIj{0) = iIj{6 ± 2k7r) and I {6) = I {6 ± 2/c7r) {k = 1, 2, 3 . . .). The total power of the probe P and its stationary mode 
with real propagation constant /3, 



(2) 



V^(i9,z) = e^^^^((9). (3) 

From the Ref. [lH, if the rotating speed a; 7^ 0, 0(6>) is a complex function. The stability of the stationary mode 
can be numerical identified by the computation of eigenvalues for small perturbations and direct simulations. The 
perturbed solution is given as iIj = e^^^[0((9) + u{e)e'^^ + v*(l9)e-^^^]. Substitution of this ansatz into Eq. ([I]) and 
linearization lead to the eigenvalue problem, 

where F = l + /+|0p. The solution is stable if all the eigenvalues of Eq. (j4]) are real. 

Firstly, we apply the voltage U along the positive direction of the extraordinary axis of BPRC, leading to £^0 > 
which results in Eq. ([1]) featuring self- focusing (SF) nonlinearity. The intensity of the induction field (alias, the 
pumps) is given as / = Asin^^. According to Eq. ([1]), the maximums of the intensity induces the minimums of the 
linear potential at ^ = ±7r/2 (The numerical domain here is selected in — 7r/2 < 9 < 7r/2), which creates the DWP 
along the circle. Therefore, by giving a set of parameter {Eq^uo^ A^P)^ we can obtained the solution from Eq. ([T]). 
The mode solutions in this work are drawn from the numerical code named "PCSOM" which is borrowed from Ref. 
[35} . For the SF case, because the symmetry breaking occurs in the GS mode, other types of mode are not considered 
here. Typical examples of the stable symmetric and asymmetric GS mode with nonzero rotating speed are displayed 
in Fig. [2j The figures imply that increasing £^0 leads to the symmetry breaking of the GS mode. The stability regions 
of these two modes drawn in the plane of (£^0,^) with different A and P in the first rotational Brillouin zone (FRBZ, 
the interv al of < u < 1/2) [3 are shown in panel |3(a)p(c)] of Fig. O 

In Fig. |3(a)||3(d)| it is interesting to find a special point (labeled as 'G' inside the panels) dividing the plane into 
three areas (symmetric area, asymmetric area and bistable area), which plays the role similar to a triple phase point. 
Prior to the point {uj < cjg), the symmetric mode transits to the asymmetric mode directly, and afterwards {u > cjg). 
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FIG. 2: (Color online) Examples of stable symmetric and asymmetric GS modes, found in the system with SF case at rotation 
speed uj = 0.25, with the parameters {Eq, A, P) = (2, 2, 0.5) and (4, 2, 0.5). Panels (a), (e) display, severally, real and imaginary 
parts of the symmetric and asymmetric GS modes, while (b), (f) show their local-power (density) profiles. The dashed-dot 
curves in (b) and (f) depict the corresponding intensity profile, 1(0); in the present case, it is sin^^. Panels (c), (g) exhibit 
the growing rates of the wave modes. And panel (d), (h) demonstrate the evolutions of these wave modes perturbed by 10% 
noises, respectively. 
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FIG. 3: Stability diagrams in the plane of (Eo^uj) inside the FRBZ. (a), (b) and (c) are for the SF case, while (d) is for the 
SDF case. In panel (a), (b) and (c), dark grey (bottom) represents the stability region of symmetric GS mode, light grey (top) 
represents the stability region of asymmetric GS mode and middle grey (right edge) represents the bistable region of these two 
types of modes. The parameters of them are (a) A = 2, P = 0.5, {h)A = 3, P = 1 and (c) A — 2, P — 0.25. In panel (d), 
dark grey (bottom) represents the stability region of antisymmetric mode, light grey (top) represents the stability region of 
antiasymmetric mode and middle grey (right edge) represents the bistable region of these two types of modes. The parameter 
of it is ^ = 3, P = 0.5. 
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FIG. 4: The ASP of the GS symmetric and asymmetric mode, defined by Eq. (|5]), as function of the Eq before (a) and after (b) 
the 'G' point. The sohd hne represents the branch where the stable solution are found, while the dot line represent the finding 
of unstable solution. Panel (a) shows the bifurcation of the supercritical type, while panel (b) shows the subcritical type. 
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FIG. 5: (Color online) (a) Coordinate uog as a function of P with A = 2,3 and 4. (b) Coordinate Eq as a function of P with 
A = 2,3 and 4. 



the transition happens after crossing the bistable area. This phenomenon implies that this point also separates the 
types of transition between the symmetric and asymmetric modes. To identified the symmetry breaking transition 
(alias SBB), we defined the asymmetry parameter (ASP) of the GS mode as 



ASP 



Jo J-TT 



\(j)\^d0 



/P. (5) 



Typical bifurcation diagrams of SBB on the either side of the 'G' point expressed by the ASP(£^o) are shown in Fig. 
m The plots demonstrating the SBB prior to and after the 'G' point are the supercritical type and subcritical type, 
respectively. The 'G' point, which has the rotating speed close to the right edge of the FRBZ, denotes the transition 
between these two kinds of SBB. The positions of point 'G' in the (£^0,^) plane are determined by the coordinates 
{Eq^ujg)- These coordinates as a function of total power P with different values of A (the intensity of the pumps) are 
displayed in Fig. O The panels in Fig. [5] show that larger power of the probe, P, leads to the transition point moving 
toward the left side of the plane; higher intensity of the pump leads to the point moving upward. These properties 
give the possibility to control the types of SBB in such rotating system. 

If we apply the voltage U along the negative direction of the extraordinary axis of BPRC, it gives rise to < 
and results in Eq. ([1]) featuring self-defocusing (SDF) nonlinearity. Under this circumstance, the maximum intensity 
of the pump becomes the maximum of the linear potential. Therefore, in the SDF case, we select the intensity profile 
of the pump as / = Acos^O^ which keeps the minimums of the linear potential appearing at = ±7r/2 (recall that 
our numerical domain is select to — tt < < it). The symmetry breakings in SDF case are occured in the non-ground 
state. In the DWP, it takes place between the antisymmetric and ant i- asymmetric state. Therefore, the symmetry 
breaking of this type is named as antisymmetry breaking. Typical examples of antisymmetric and anti-asymmetric 
modes are displayed in Fig. [6l However, the stability plane of (Eq^uj) in panel [3(d)] shows no transition point is found. 
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FIG. 6: (Color online) Examples of stable antisymmetric and anti-asymmetric modes, found in the system with SDF case at 
rotation speed uj = 0.25, with the parameters {Eo,A,P) = (2,2,0.5) and (4,2,0.5). Panels (a), (c) display, severally, real and 
imaginary parts of the antisymmetric and anti-asymmetric modes, while (b), (d) show their local-power (density) profiles. The 
dashed-dot curves in (b) and (d) depict the corresponding intensity profile, 1(0); in the present case, it is cos^ 0. Panels (c), (g) 



exhibit the growing rates of the wave modes, 
by 10% noises, respectively. 



And panel (d), (h) demonstrate the evolutions of these wave modes perturbed 



The antisymmetric mode transit to the ant i- asymmetric mode directly, implying the type of SBB of antisymmetry 
breaking is only supercritical @. Therefore, in the SDF case, we cannot switch the type of SBB. 

Conclusion: the objective of this work is to study the possibility of switching (or controlling) the type of SBB 
in the rotating DWP system made onto the BPRCs. The rotating potential can be created by means of optical- 
induction techniques, using vortex pump waves with the ordinary polarization in a photorefractive material. Increase 
the magnitude of the biased voltage can lead to the symmetry breaking of the mode inside the waveguide. The 
numerical analysis of the symmetry breaking is studied inside the FRBZ. For SF case, the symmetry breaking occurs 
in the GS. A special triple-like point, which separates the type of SBB (supercritical and subcritical), is found close 
to the right edge of the FRBZ. On the left side of the point, the type of SBB is supercritical, while on the right side 
of the point it is subcritical. The influence on the position of the point in the (^o,^) plane is also studied. It shows 
that the type of SBB can be switched in the rotating system in the SF case. For the SDF case, the type of symmetry 
breaking is antisymmetry breaking, which takes place between the antisymmetric and anti-asymmetric mode. The 
type of SBB of the antisymmetry breaking is only the supercritical type, which results in the switch of SBB cannot 
be achieved in this case. 

This work was support by the CNNSF (grant No. 11204089,11205063). 



[1] L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Nauka Publishers, Moscow, 1974). 

[2] P. G. Kevrekidis, Z. Chen, B. A. Malomed, D. J. Frantzeskakis, and M. I. Weinstein, Phys. Lett. A 340 (2005) 275. 
[3] G. Theocharis, P. G. Kevrekidis, D. J. Frantzeskakis, and P. Schmelcher, Phys. Rev. E 74 (2006) 056608. 
[4] L. Albuch, B. A. Malomed, Math. Comput. Simulat. 74 (2007) 312. 

[5] G. Herring, P. G. Kevrekidis, B. A. Malomed, R. Carretero- Gonzalez, D. J. Frantzeskakis, Phys. Rev. E 76 (2007) 066606. 
[6] M. Matuszewski, B. A. Malomed, and M. Trippenbach, Phys. Rev. A 75 (2007) 063621. 

[7] M. Trippenbach, E. Infeld, J. Gocalek, M. Matuszewski, M. Oberthaler, and B. A. Malomed, Phys. Rev. A 78 (2008) 
013603. 

[8] I. I. Satija, R. Balakrishnan, P. Naudus, J. Reward, M. Edwards, and C. W. Clark, Phys. Rev. A 79 (2009) 033616. 

[9] B. Xiong, J. Gong, H. Pu, W. Bao, and B. Li, Phys. Rev. A 79 (2009) 013626. 
[10] T. Mayteevarunyoo, B. A. Malomed, and G. Dong, Phys. Rev. A 78 (2008) 053601. 
[11] L. Salasnich, B. A. Malomed, and F. Toigo, Phys. Rev. A 81 (2010) 045603. 
[12] S. K. Adhikari, B. A. Malomed, L. Salasnich, and F. Toigo, Phys. Rev. A 81 (2010) 053630. 
[13] N. Bar-Gill, G. Kurizki, B. A. Malomed, and C. Sudheesh, Phys. Rev. A 82 (2010) 013610. 



6 



[14] H. Sakaguchi and B. A. Malomed,Phys. Rev. E 83 (2011) 036608. 

[15] V. A. Brazhnyi, and B. A. Malomed, Phys. Rev. A 83 (2011) 053844. 

[16] A. R. Davoyan, I. V. Shadrivov, and Y. S. Kivshar, Opt. Lett. 36 (2011) 930. 

[17] Y. Li, W. Pang, S. Fu, and B. A. Malomed, Phsy. Rev. A 85, (2012) 053821. 

[18] Y. Li, W. Pang, B. A. Malomed, Phys. Rev. A 86 (2012) 023832; e-print larXiv: 1209.3166 1 (2012) . 

[19] G. Kakarantzas, A. Ortigosa- Blanch, T.A.Birks, P. St. J. Russell, L. Farr, F. Couny, and B. J. Mangan, Opt. Lett. 28 

(2003) 158. 

[20] H. Saito and M. Ueda, Phys. Rev. Lett. 93 (2004) 220402. 

[21] S. Schwartz, M. Cozzini, C. Menotti, L Camsotto, P. Bouyer, and S. Stringari, New J. Phys. 8 (2006) 162. 
[22] Y. V. Kartashov, B. A. Malomed, and L. Torner, Phys. Rev. A 75 (2007) 061602(R). 
[23] S. Jia and J. W. Fleischer, Phys. Rev. A 79 (2009) 041804(R). 
[24] A. E. Fetter, Rev. Mod. Phys. 81 (2009) 647. 

[25] L. Wen and X. Luo, Laser Phys. Lett. 9 (2012) 618; L. Wen, H. Xiong, and B. Wu, Phys. Rev. A 82 (2010) 053627. 
[26] Hidetsugu Sakaguchi and Kenji Kita, J. Phys. Soc. Jpn. 81 (2012) 104008. 

[27] N. K. Efremidis, S. Sears, and D. N. Christodoulides, J. W. Fleischer and M. Segev, Phys. Rev. E 66 (2002) 046602; J. 

W. Fleischer, M. Segev, N. K. Efremidis, D. N. Christodoulides, Nature (London) 422 (2003) 147. 
[28] J. W. Fleischer, G. Bartal, O. Cohen, O. Manela, M. Segev, J. Hudock, and D. N. Christodoulides, Phys. Rev. Lett. 92 

(2004) 123904. 

[29] H. Buljan, O. Cohen, J. W. Fleischer, T. Schwartz, and M. Segev, Z. H. Musslimani, N. K. Efremidis, and D. N. 

Christodoulides, Phys. Rev. Lett. 92 (2004) 223901; O. Cohen, G. Bartal, H. Buljan, T. Carmon, J. W. Fleischer, M. 

Segev, D. N. Christodoulides, Nature (London), 433 (2005) 500. 
[30] X. Wang, Z. Chen and P. G. Kevrekidis , Phys. Rev. Lett. 96 (2006) 083904; X. Wang, Z. Chen, J. Yang, 31 (2006) 1887. 
[31] F. Lederer, G. L Stegemanb, D. N. Christodoulides, G. Assanto, M. Segev, Y. Silberberg, Phys. Rep. 463 (2008) 1. 
[32] P. Zhang, S. Liu, J. Zhao, C. Lou, J. Xu, and Z. Chen, Opt. Lett. 33 (2008) 878. 
[33] Y. Li, W. Pang, Y. Chen, Z. Yu, J. Zhou, and H. Zhang, Phys. Rev. A 80 (2009) 043824. 
[34] Z. Chen, M. Segev and D. N Christodoulides, Rep. Prog. Phys. 75 (2012) 086401. 
[35] J. Yang and T. L Lakoba, Stud. Appl. Math. 118 (2007) 153; 120 (2008) 265. 



